Spin Hall effect of alloys: Extrinsic and intrinsic contribution 
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A fully relativistic description of the spin-orbit induced spin Hall effect is presented that is based 
on Kubo's linear response formalism. Using an appropriate operator for the spin current density 
a Kubo-Stfeda-like equation for the spin Hall conductivity (SHC) is obtained. An implementa- 
tion using the Korringa-Kohn-Rostoker (KKR) band structure method in combination with the 
Coherent Potential Approximation (CPA) allow detailed investigations on various alloy systems. A 
decomposition of the SHC into intrinsic and extrinsic contributions is suggested. Accompanying 
calculations for the skew-scattering contribution of the SHC using the Boltzmann equation demon- 
strate the equivalence to the Kubo formalism in the dilute alloy regime and support the suggested 
decomposition scheme. 
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The emerging research field spintronics has been de- 
veloped very rapidly during the last years. The reason 
for the broad interest in this field is based on the close 
connection to fundamental scientific questions as well as 
its impact on technology [J, 12]. In this context, the spin 
Hall effect (SHE) is one of the most promising phenom- 
ena. It denotes the observation that a charge current 
flowing through a solid is accompanied by a transversal 
spin current. This occurs even for non- magnetic solids as 
was demonstrated by experiments on pure Pt [3j]. 

Both the anomalous Hall effect (AHE) in ferromagncts 
and the SHE are caused by the influence of spin-orbit cou- 
pling (SOC). Accordingly, their theoretical description is 
quite similar 0-0 • For ideal systems an intrinsic 
mechanism was identified that allows to express the cor- 
responding response function in terms of the Berry cur- 
vature [f^ HI . On this basis ah initio calculations for the 
intrinsic spin Hall conductivity (SHC) were performed 
[jaWl^. As for the AHE, the additional extrinsic SHC in 
dilute and concentrated alloys is ascribed to skew and 
side-jump scattering caused by SOC. The role of these 
mechanisms for the SHE has been studied so far primar- 
ily by model calculations [3, E3] • First principle calcula- 
tions for the extrinsic SHC of dilute alloys on the basis of 
the Boltzmann formalism that account for the skew scat- 
tering mechanism have been performed only very recently 

15l[l6|. However, a complete description of intrinsic and 
extrinsic mechanisms giving rise to the SHE applicable 
to ideal as well as alloy systems, as it is presented below, 
was missing so far. As pointed out by several authors 

17 . HH, a central issue for such an approach is an ad- 
equate definition for the spin current density operator 
that accounts for SOC. This was supplied recently by 
Vernes et al. [l9| by starting from the Bargmann-Wigner 
four- vector spin polarization operator T [20]. Demand- 
ing that the spin polarization is connected with the spin 



current density via a corresponding continuity equation 
an explicit expression for the spin current operator was 
given. 

An adequate formal basis for the discussion of the SHE 
in non-magnetic metals is supplied by Kubo's linear re- 
sponse formalism that allows to derive an expression for 
the spin Hall conductivity tensor. To avoid any approxi- 
mation when dealing with SOC the underlying electronic 
structure is described in a fully relativistic way by the 
four-component Dirac formalism 21]. As for the ordinary 
electrical conductivity the perturbation due to the exter- 
nal electric field is represented in terms of the current 
density operator j. In its relativistic form this operator 
is given by: 



e col 
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where a is the vector-matrix of the standard Dirac ma- 
trices at 2l| and the other quantities have their usual 



meaning. The response function to be considered for the 
SHE is the spin current density. Considering for the z- 
component of the spin polarization vector the current 
density along the x direction the corresponding operator 
is given by [lj| El ^ : 



eca x 
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where p z is the canonical momentum operator, E z is the 
z component of the vector of the relativistic spin matrix, 
(3 and 75 are Dirac matrices [211 ] . 

Adopting a single-particle description of the electronic 
structure in terms of the retarded (G + ) and advanced 
(G~) Green's function and restricting to T = K, an 
explicit expression for the SHC is obtained that is simi- 
lar to the Kubo-Stfeda equation for the anomalous Hall 
conductivity cr xy of ferromagnetic systems. Considering 
for the spin polarization along z its current density along 



2 



x due to an electric field along y, the SHC a\ y * s given 
bydl: 

h 



crL, = 



Xy 27TiVfi 



Tr <J|G+j y G-) c 



where terms containing products of the retarded (or ad- 
vanced) Green's functions have been dropped [6\. Due 
to symmetry the last term is site-diagonal for the cubic 
systems considered here. As furthermore all system con- 
sidered here are metallic 0, EH it has been omitted as 
well. 

The electronic Green's function occurring in Eq. (|3]) 
for the Fermi energy Ep can be evaluated in a very effi- 
cient way by use of the relativistic version of the multiple 
scattering or Korringa-Kohn-Rostoker (KKR) formalism 
pBf . This approach is applicable to ideal systems but 
can also be applied straightforwardly to disordered al- 
loys. For this case the brackets (...) c in Eq. <j3j> imply 
a configurational average that is performed within the 
Coherent Potential Approximation (CPA) [26j]. For this 
purpose the KKR-CPA approach used for the electrical 
conductivity tensor of alloys on the basis of the Kubo- 
Stfeda equation 27J has been adapted to Eq. Q. In 
particular, this approach accounts explicitly for the so- 
called vertex corrections those represent the difference in 
the correlated and uncorrelated configurational averages 
of the type {J z G + j y G~) c and (J^G + ) c (j y G^) c , respec- 
tively. For the following it is important to note that the 
vertex corrections correspond to the scattering-in term 
within the Boltzmann formalism [26j j . 

Representing the anomalous Hall conductivity (AHC) 
<7 xy in terms of Feynman diagrams it was demonstrated 
that all extrinsic contributions to a xy due to skew and 
side-jump scattering correspond to terms involving ver- 
tex corrections [281 ] - Obviously, the same conclusion can 
be drawn for the SHC crj y . The remaining diagrams, 
involving no vertex corrections, stand for the intrinsic 
anomalous or spin Hall conductivity, plus corrections to 
this due to chemical disorder. It therefore seems natu- 
ral to extend the definition of the intrinsic SHC a 
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to 



diluted and concentrated alloys to represent all contribu- 
tions not connected to the vertex-corrections. According 
to this definition, Eq. ©, including the vertex correc- 
tions, gives the total SHC a\ y = a xy C while the intrinsic 
SHC a^ y ntI = er^y oVC is obtained if those are ignored. 
Thus, a z J' ltT is the intrinsic SHC of the effective CPA 
medium specific for each composition of a certain alloy. 
As a consequence, the extrinsic SHC a^. y xt1 ' to be ascribed 
to the skew and side-jump scattering mechanisms is ob- 
tained from the difference a^ y xtr = a* y — cr* y ntr . 

The approach sketched above has been applied to in- 
vestigate the SHE for the fee alloy systems Au x Pti_a; 
and Ag^Aui-a;. As can be seen from Fig. [T] cr^ y ntr ob- 
tained from Eq. ([3]) ignoring the vertex corrections varies 
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FIG. 1: (Color online) Spin Hall conductivity o\ y for the 
alloy systems Au^Pti-^ and Ag x Au\- x . The full squares 
correspond to calculations including vertex corrections while 
the open squares represent calculations without vertex cor- 
rections. Results from other ab initio investigations on the 
intrinsic SHE of pure Pt HH and Au HQ2 

are included. Be- 
cause of the different definitions for the spin current operator 
the results from Refs. [Tl] and have been multiplied by a 
factor of 2 for the sake of consistency. 



nearly linearly with the concentration throughout the 
whole composition regime. In addition, Fig. H] (top) 
shows results for the intrinsic SHC of Au [9|, [l2| and 
Pt [11] obtained by other authors using an expression 
for (Txy ntr i n terms of the Berry curvature. Taking into 
account the differences between the various calculation 
schemes used - in particular concerning the treatment of 
spin-orbit coupling and the definition of the spin current 
density operator - these data fit reasonably well with 
the alloy data obtained using the Kubo-Stfeda-like equa- 
tion (Eq. (|3])). This obviously justifies the extension of 
the definition for o-^ y ntr to the alloy case to represent all 
terms that do not involve the vertex corrections. 

In contrast to the intrinsic SHC, the total one (cr^y) 
shows for both investigated alloy systems a divergent 
behavior in the dilute limit when the concentration x 
approaches or 1, respectively. Interestingly, the cor- 



responding extrinsic SHC cr%. y * tr changes sign when the 
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concentration x varies from to 1. For concentrated al- 
loys (0.2 < x < 0.8) the intrinsic and total SHC do not 
differ strongly. As this behavior is also found for other 
alloy systems that do not show a change in sign for the ex- 
trinsic SHC er^y Xtr it seems that the impact of the vertex 
corrections in the concentrated alloy regime is in general 
negligible. 

For the anomalous Hall effect it is known that one may 
classify solid state materials according to the scaling rela- 
tion between the AHC cr xy and the longitudinal conduc- 
tivity <t xx . For metallic systems those fall in general into 
the so-called ultra-clean regime (cr xx > 10 6 (ftem) -1 ) the 
skew-scattering mechanism should dominate a 



In this case the scaling relation cr xy « cr X y GW = S cr xx holds 
with S the so-called skewness factor. Assuming the same 
to apply for the SHE as well, one may expect for the ex- 
trinsic SHC the relation: 



z cxtr z skew 

x y x y 
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where S z is the corresponding skewness factor and the 
term <r X y J represents the contribution due to the side- 
jump mechanism. 

Plotting the extrinsic SHC cr X y Xtr of AuzPti-z and 
Ag^Aui-z versus the corresponding er xx with the concen- 
tration as an implicit parameter indeed a linear behavior 
is found in the dilute regimes (x < 0.1 or x > 0.9), as can 
be seen from Fig. [21 Fitting a straight line to the data 
for the considered systems and extrapolating to er xx = 
allows to deduce the side-jump contribution on the ba- 
sis of Eq. ([4|. The results obtained for Au a; Pti_ a ; and 
Ag^Aui-a; are shown in Fig. [3] together with data ob- 
tained for two other alloy systems. Obviously, may 
take either sign and is in most cases found to be much 
smaller than the intrinsic contribution. On the other 
hand, a Z y kow is dominating in the dilute limit since it 
scales with the concentration, while cr^ y ntr does not de- 
pend on the concentration explicitly. However, an actual 
impurity concentration for the crossover between intrin- 
sic and skew scattering regime depends on the considered 
alloy. 

To support the analysis of the results for the total 
SHC presented above, complementary work has been 
done using the Boltzmann formalism for the SHE fljj ]. 
The results for er xx obtained this way are found in very 
good agreement with those obtained using the Kubo- 
Greenwood equation (see top panel of Fig. |4|). As men- 
tioned above, the vertex corrections giving rise to <7 X y Xtr 
correspond to the scattering-in processes occurring in 
the Boltzmann formalism. As it was demonstrated re- 
cently, the latter ones give rise to the skew-scattering 
mechanism [l^|. Corresponding results for <7 X y kew (bot- 
tom panel of Fig. (jU) are also found in very satisfying 
agreement with the results based on the Kubo-Stfeda-like 
formula (Eq. (|3J|) together with the described decompo- 
sition. This finding convincingly shows the equivalence 
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FIG. 2: (Color online) The extrinsic spin Hall conductivity 
°xy X 1 versus the longitudinal conductivity cr xx for Ag^Aui-x 
and Au^Pti-! (black line/circles). The blue and orange lines 
are explained in the text. 
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FIG. 3: (Color online) The side-jump and intrinsic contri- 
butions to the spin Hall conductivity, <r X y J and o- X y ntr , respec- 
tively, for various dilute transition metal alloys A(B) with the 
concentration of 1 at.% for the dissolved component B. 
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"Spinphanomene in reduzierten Dimensionen" . Further 
more, this work was supported by the International Max 
Planck Research School for Science and Technology. 
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FIG. 4: (Color online) The longitudinal conductivity <r xx 
(top) and the skew scattering contribution to the spin Hall 
conductivity <r X y kcw (bottom) for various dilute alloys A(B) 
obtained on the basis of the Kubo-Stfeda equation (left bar) 
and the Boltzmann formalism (right bar), respectively. The 
concentration for the dissolved component B is 1 at.%. 



of both approaches for the dilute alloy regime and it also 
justifies once more the used definition for the intrinsic 
SHC introduced above. 

In summary, an expression for the spin Hall conductiv- 
ity <7 xy has been derived in analogy to the Kubo-Stfeda 
equation for the anomalous Hall conductivity cr xy of fer- 
romagnets. An implementation within the fully relativis- 
tic KKR-CPA formalism allows numerical investigations 
for various transition metal alloy systems over the whole 
range of composition. We decompose the total SHC a^ y 
into its intrinsic and extrinsic parts. In the concentrated 
alloys the intrinsic contribution of the effective medium 
always dominates. The extrinsic contribution, on the 
other hand, shows in general a diverging behavior in the 
dilute alloy regime that is ascribed to the skew scatter- 
ing contribution. Accompanying calculations on the basis 
of the Boltzmann formalism demonstrate its equivalence 
with the Kubo formalism in the dilute alloy regime and 
support the decomposition made for the total SHC u xy . 

The authors S.L., D.K. and H.E. would like to thank 
the DFG for financial support within the SFB 689 
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